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Much effort in cancer research has focused on the tiny part of our genome that codes for mRNA. However, it has recently 
been recognized that microRNAs also contribute decisively to tumorigenesis. Studies have also shown that epigenetic 
silencing by CpG island hypermethylation of microRNAs with tumor suppressor activities is a common feature of human 
cancer. The importance of other classes of non-coding RNAs, such as long intergenicncRNAs (lincRNAs) and transcribed 
ultraconserved regions (T-UCRs) as altered elements in neoplasia, is also gaining recognition. Thus, we wondered 
whether there were other ncRNAs undergoing CpG island hypermethylation-associated inactivation in cancer cells. We 
focused on the small nucleolar RNAs (snoRNAs), a subset of ncRNA with a wide variety of cellular functions, such as 
chemical modification of RNA, pre-RNA processing and control of alternative splicing. By data mining snoRNA databases 
and the scientific literature, we selected 49 snoRNAs that had a CpG island within < 2 Kb or that were processed from 
a host gene with a 5'-CpG island. Bisulfite genomic sequencing of multiple clones in normal colon mucosa and the 
colorectal cancer cell line HCT-116 showed that 46 snoRNAs were equally methylated in the two samples: completely 
unmethylated (n = 26) or fully methylated (n = 20). Most interestingly, the host gene-associated 5'-CpG islands of the 
snoRNAs SNORD123, U70C and ACA59B were hypermethylated in the cancer cells but not in the corresponding normal 
tissue. CpG island hypermethylation was associated with the transcriptional silencing of the respective snoRNAs. Results 
of a DNA methylation microarray platform in a comprehensive collection of normal tissues, cancer cell lines and primary 
malignancies demonstrated that the observed hypermethylation of snoRNAs was a common feature of various tumor 
types, particularly in leukemias. Overall, our findings indicate the existence of a new subclass of ncRNAs, snoRNAs, that 
are targeted by epigenetic inactivation in human cancer. 



Introduction 

Coding exons account for only 1.5% of the genome, 1 despite 
being the focus of most of the current biomedical research. A 
large proportion of the genome is made up of non-protein cod- 
ing regions that might have critical biological important, contain- 
ing gene regulatory regions (transcriptional and splicing types), 
matrix attachment sites, origins of replication, other functional 
elements and non-coding RNAs (ncRNAs). 2,3 The physiological 
and pathological importance of this functional part of the non- 
protein-coding genome is particularly apparent in a large class 
of small non-coding RNAs (ncRNAs) known as microRNAs. 4 
These are about 22 nucleotides long and repress gene expression 
in a variety of eukaryotic organisms. 4 In human cancer, miRNA 
expression profiles differ between normal tissues and derived 
tumors and between tumor types. 5 miRNAs can act as onco- 
genes or tumor suppressors, with a key role in tumorigenesis. 6 " 8 
Defects in miRNA function have been associated with a failure of 



Correspondence to: Manel Esteller; Email: mesteller@idibell.cat 
Submitted: 12/05/11; Revised: 01/11/12; Accepted: 01/13/12 
http://dx.doi.Org/1 0.41 61/rna.l 9353 






www.landesbioscience.com 


RNA Biology 


881 



miRNA post-transcriptional regulation, 5 miRNA transcriptional 
repression by oncogenic factors, 10 loss-of-function genetic altera- 
tions in the genes involved in miRNA-processing pathways"" 13 
and transcriptional silencing associated with hypermethylation 
of CpG island promoters. 14 Thus, as occurs with miRNAs, it is 
likely that other types of ncRNAs are also involved in human 
tumorigenesis and are characterized by epigenetic and genetic 
defects in this disease. 3 In this context, Ultraconserved Regions 
(UCRs), a subset of conserved sequences that are located in intra- 
genic and intergenic regions, 15,16 are altered at the transcriptional 
level in human tumorigenesis. 17 Interestingly, transcribed UCRs 
(T-UCRs) undergo DNA methylation-associated silencing in 
cancer cells. 18 

Another important class of ncRNAs that are potentially 
altered in human cancer are the small nucleolar RNAs (snoR- 
NAs), which are localized in the nucleolus. 19 They are responsible 
for methylation 20,21 and pseudouridylation 22,23 of rRNA (rRNA) 
at about 50-100 sites per eukaryotic ribosome. However there 



is increasing evidence of the targeting of other classes of RNAs, 
such as mRNAs. 24 snoRNAs are divided into two main classes: 
box C/D and box H/ACA, 25 on the basis of their conserved sec- 
ondary structural characteristics and associated modification 
reactions. 19,26 The C/D box snoRNAs guide the position-specific 
2'-0-methylation and are associated with a core of four proteins: 
fibrillarin (the methyltransferase), NOP56, NOP5/NOP58 and 
NHPX. The H/ACA snoRNAs direct RNA pseudouridylation 
of rRNA and are associated with dyskerin (the pseudouridine 
synthase), GAR1, NHP2 and NOP10. 24,26,27 Mutations in the 
human dyskerin gene, NHP215 and NOP1016 gene are associ- 
ated with the X-linked genetic disorder, dyskeratosiscongenita 
(DC), where malfunction of rRNA and shortening of telomeres 
have been observed. 28 " 30 As mutations in the dyskerin gene have 
also been associated with cancer susceptibility, 28 " 30 it was sug- 
gested that snoRNAs were involved in the onset and progres- 
sion of cancer. One of the first studies to address this possibility 
reported a snoRNA to be highly expressed in normal brain, but 
significantly reduced in meningioma. 31 Other studies showed 
that GAS5, a snoRNA-host gene, controls cell survival by induc- 
ing or sensitizing cells to apoptosis. 32,33 A substantial decrease of 
GAS5 in breast cancer samples compared with adjacent unaf- 
fected normal breast epithelial tissues also suggests that it has a 
role as a tumor suppressor gene. 33 

The association between snoRNAs and cancer was underlined 
by other studies showing that a homozygous 2-bp (TT) deletion 
of the snoRNA U50 is strongly associated with prostate cancer 34 
and that U50 undergoes frequent genomic heterozygous deletion 
and transcriptional downregulation in breast cancer. 35 U50 over- 
expression reduces colony-forming ability in prostate and breast 
cancer cells. 34,35 Taken together, these studies suggest that non- 
coding snoRNA U50 is important in the development and pro- 
gression of breast and prostate cancers. 34,35 More recently, it was 
shown that a diversity of snoRNAs are differentially expressed 
in non-small-cell lung cancer with respect to the corresponding 
matched tissue, 36 encouraging investigation into the possible role 
of snoRNAs in oncogenesis. Another study has linked at least 
one snoRNA to the post-transcriptional processing of a protein- 
coding gene. 37 There is also evidence that other snoRNAs might 
be involved in the regulation of gene expression by giving rise to 
other regulatory RNA species, such as miRNAs, linking snoR- 
NAs to RNA silencing. 38 It would therefore be very interesting to 
identify the function and mechanisms of the orphan snoRNAs. 

The downregulation of tumor suppressor protein-coding 
genes (e.g., hMLHl, BRCA1, VHL and pl6 INK4a ), 39,40 and 
ncRNAs with growth inhibitory functions, such as miRNAs 14 
has been closely linked to the presence of CpG island promoter 
hypermethylation. Thus, we wondered whether the same mecha- 
nism could play a role in the loss of adequate snoRNA expres- 
sion in tumors. Usually, snoRNAs are genomically found in the 
introns of protein-coding or non-protein-coding host genes, with 
each intron carrying only one single snoRNA and their transcrip- 
tion being synchronized with that of the host gene. After splicing 
they are generally processed by debranching and exonucleolytic 
trimming of the 5'- and 3'-ends, 24,41 " 43 and assembled with spe- 
cific core proteins that are essential for their localization and 



correct enzymatic activity, and for preventing their degradation. 27 
However, intergenicsnoRNAs are independently transcribed by 
RNA polll as independent units, 24 and some human intron- 
encoded snoRNAs may have their own independent promoters. 23 
Herein, we present a double candidate and genomic approach to 
unmask snoRNA-associated CpG islands that undergo cancer- 
specific hypermethylation-associated transcriptional silencing, 
such as SNORD123, U70C and ACA59B. These findings sup- 
port a model in which epigenetic disruption of emerging new 
classes of ncRNAs, such as snoRNAs, is a common feature of 
human tumorigenesis. 

Results and Discussion 

snoRNA CpGisland DNA methylation analyses. To identify 
snoRNAs with putative DNA methylation-related inactivation in 
human tumors, we data-mined the scientific literature on snoR- 
NAs published during 2000-2011, as made available by PubMed. 
gov, the human genome browser at UCSC 44 and the snoRNA- 
LBME-db database. 24 The CpG Island Searcher Program 45 was 
used to determine which snoRNAs were located within <2 Kb of 
a CpG island, since it has been estimated that more than 90% 
of the human promoters of another type of ncRNA, the microR- 
NAs, are located 1 Kb upstream of the mature transcript. 46,47 
The DNA methylation status of CpG islands within 2 Kb are 
also important for regulating the expression of a second type of 
ncRNA, the T-UCRs. 18 We also included snoRNAs that were 
processed from a host gene RNA containing a 5 -CpG island. 
Figure 1 illustrates both categories of snoRNA-related CpG 
islands. Using the described conditions, we selected 49 snoRNAs 
that had a CpG island within a distance of < 2 Kb (15 intergenic 
independent snoRNAs and 24 within an intron of a host gene) 
or that were processed from the transcript of a host gene with a 
5'-CpG island (n = 10). The characteristics of the 49 selected 
snoRNAs and the summarized results are shown in Table 1. 

We performed bisulfite genomic sequencing of multiple clones 
using primers encompassing the 49 snoRNA-associated CpG 
islands to determine the DNA methylation patterns in normal 
colon mucosa and the colorectal cancer cell line HCT-116. We 
observed a completely unmethylated status for 23 snoRNA- 
related CpG islands in normal tissue and colon cancer cells 
(Table 1). Examples of the bisulfite sequencing analyses are 
shown in Figure 2. We also found a dense DNA methylation 
profile for 20 snoRNA-associated CpG islands in normal colon 
mucosa and HCT-116 colorectal cancer cells (Table 1). Examples 
of the bisulfite sequencing analyses are shown in Figure 3. Most 
importantly, we found a cancer-specific hypermethylation event 
for the snoRNAs SNORD123, U70C and ACA59B. Their asso- 
ciated CpG islands were completely unmethylated in normal 
colon mucosa and heavily hypermethylated in HCT-116 colorec- 
tal cancer cells (Table 1 and Fig. 4). For all three cases, the CpG 
islands studied were in the 5'-region of the host gene where the 
snoRNA is located: the long non-coding gene LOC100505806 
(SNORD123), Astrotactin 2 (U70C) and Solute Carrier Family 
47 Member 1 (ACA59B). The DNA methylation results were also 
confirmed using methylation-specific PCR (Fig. SI). 
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Hypermethylation of snoRNA-related CpG islands is associ- 
ated with transcriptional silencing. To demonstrate transcrip- 
tional silencing of these snoRNAs in cancer cells in association 
with the presence of CpG island hypermethylation, we measured 
transcript levels by quantitative RT-PCR. For the SNORD123, 
we analyzed the expression of the snoRNA itself, the long non- 
coding RNA LOC100505806 from which the snoRNA is pro- 
cessed and the mRNA of the Semaphorin 5A (SEMA5A) that it 
is transcribed in the opposite direction from the same CpG island 
(Fig. 4). No expression of the SNORD123, LOC100505806 
and SEMA5A transcripts could be detected in HCT-116 cells 
in which the shared CpG island was methylated (Fig. 5A). 
Normal colon mucosa expressed the three transcripts (Fig. 5A). 
Methylation-specific PCR analyses of two additional colon can- 
cer cell lines identified a hypermethylated and unmethylatedCpG 
island in SW48 and DLD1 cells, respectively (Fig. SI). Loss of 
the SNORD123, LOC100505806 and SEMA5A transcripts was 
observed in the hypermethylated SW48 cells and expression of 
the three transcripts was found in the unmethylated DLD1 cells 
(Fig. 5A). The absence of the SNORD123 transcript in HCT-116 
and SW48 cells and its presence in DLD1 cells was also validated 
by northern-blot analyses (Fig. 5B). For the snoRNAs U70C 
and ACA59B, no expression of the snoRNA U70C, its host 
gene ASTN2 and the host gene of ACA59B (SLC47A1) could 
be detected in HCT-116 cells in which the corresponding CpG 
island was methylated (Fig. S2). Normal colon mucosa expressed 
the three transcripts (Fig. S2) . Most importantly, the expression 
for ACA59B (SLC47A1) was restored upon treatment with the 
DNA demethylating agent 5-aza-2'-deoxycytidine in the HCT- 
116 cell line (Fig. S2). These results were confirmed using an 
alternative model of an isogenic HCT-116 cell line in which the 
two major DNA methyltransferases, DNMTl and DNMT3b, 
had been genetically disrupted (HCT116 DKO). 48 The CpG 
island for ACA59B (SLC47A1) was hypomethylated in DKO 
cells (Fig. S2), but hypermethylated in the HCT116 parental cell 
line. ACA59B (SLC47A1) expression was restored in DKO cells 
(Fig. S2), reinforcing the link between CpG island hypermethyl- 
ation and snoRNA silencing. 

Profile of snoRNA hypermethylation in different tumor 
types. To address the cancer-specific hypermethylation event in 
snoRNAs, we ruled out the possible presence of SNORD123, 
U70C and ACA59B CpG island methylation in a panel of nor- 
mal tissues from the colon (n = 5), breast (n = 7) and lung (n = 
22), in addition to four blood samples from healthy volunteers, 
using a DNA methylation microarray approach (Fig. 6). 49 The 
observed high-throughput DNA methylation platform included 
five, two and four CpG sites corresponding to the bisulfite 
genomic sequenced CpG islands for SNORD123, U70C and 
ACA59B, respectively (Fig. 6). The presence of SNORD123, 
U70C and ACA59B cancer-specific CpG island hypermethyl- 
ation and transcriptional silencing was not a unique feature of the 
colorectal cancer cell line HCT-116; analyzing the NCI60 panel 
of human cancer cell lines (n = 60) from nine tumor types, we 
also found them in other colon cancer cell line and lung, breast, 
prostate, ovarian, renal, melanoma, lymphoma and leukemia 
cells (Fig. 7). 
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Figure 1. Types of CpG islands associated with snoRNAs in this study. 
(A) Upstream CpG islands of a snoRNA located within 2 Kb. It includes 
intergenic independent snoRNAs and snoRNAs within an intron of a 
host gene. (B) 5'-CpG islands of host genes for which RNA processing 
generates the expression of an intronic resident snoRNA. 



Most importantly, the CpG island hypermethylation of 
SNORD123, U70C and ACA59B was not an in vitro phe- 
nomenon. Having noted the high frequency of CpG island 
hypermethylation for the three snoRNAs in leukemia cell lines 
(Fig. 7), we examined the presence of these epigenetic events in 
48 primary samples from acute lymphoblastoid leukemia (ALL) 
patients, 43 had a B-cell phenotype while five were T-ALLs 
(Fig. 8). Using the described DNA methylation platform, we 
observed SNORD123, U70C and ACA59B hypermethylation 
in 27% (13 of 48), 39% (19 of 48) and 29% (14 of 48) of acute 
lymphoblastoid leukemia cases, respectively (Fig. 8). We did not 
observe any association between CpG Island hypermethylation 
of the three studied snoRNAs and disease-free survival or overall 
survival in these patients (Log rank Mantel-Cox test p > 0.05). 
We also extended the analyses of snoRNA DNA methylation to 
primary acute myelogenous leukemia (AML) samples. Among 
16 primary acute myelogenous leukemia cases, we observed 
SNORD123 CpG island hypermethylation in 25% (4 of 16) of 
samples, while ACA59B and U70C were unmethylated in all 
cases. Finally, among 20 primary multiple myeloma cases, we 
observed ACA59B CpG island hypermethylation in 10% (2 of 
20) of samples, while SNORD123 and U70C were unmethylated 
in all cases. 

Overall, our results reveal the existence of cancer-specific 
hypermethylation events in CpG islands associated with snoR- 
NAs that lead to their transcriptional inactivation in transformed 
cells. Despite our increasing knowledge about the biological 
roles of snoRNAs, one of the main challenges in cancer research 
into ncRNAs is the identification of a particular function that 
is relevant for cellular transformation. As coding genes, 39,40 
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Table 1. DNA methylation profile of CpG islands associated with snoRNAs 



Name 


Class 


Chr. 


Strand 


Length (bp) 


Target RNA 


Host gene 


Analyzed CpG Island 


Normal Colon 


HC T116 


ACA24 


HAcaBox 


4 


+ 


131 


1SS rRNA U863 and 13S rRNA U609 


SNHG8 


CpG island within £ 2Kb 






U50 


CDBox 


G 




74 


23S rRNA C2348 and 28S rRNA G28S3 


SNHG5 


CpG island within £ 2Kb 






U5DB 


CDBox 


6 




71 


Unknown 


SNHG5 


CpG island within £ 2Kb 






ACA9 


HAcaBox 


7 




133 


28S rRNA U1670 and 23S rRNA U17S9 


C7orf40 


CpG island within £ 2Kb 






HBII-33G 


CDBox 


7 


* 


74 


18S rRNA A576 


Independent 


CpG island within £ 2Kb 






hTR 


scaRna 


3 




648 


Unknown 


Independent 


CpG island within £ 2Kb 






U91 


scaRna 


1 3 




33 


U4 snRNA C8 


Independent 


CpG island within £ 2Kb 






mgU12-22/U4-E 


scaRna 


18 


* 


421 


U4 snRNA C8 and U12 snRNA G22 


Independent 


CpG island within £ 2Kb 






HBII-382 


scaRna 


1 


+ 


82 


U2 snRNA C61 and U2 snRNA G11 


Independent 


CpG island within £ 2Kb 






mgU2-25/61 


scaRna 


1 




4/0 


U2 snRNA G25 and U2 snRNA C61 


Independent 


CpG island within £ 2Kb 






U13 


CDBox 


8 


* 


104 


Unknown 


Independent 


CpG island within £ 2Kb 






ACA62 


HAcaBox 


17 


* 


1 J3 


18S rRNA U 34 and 18S rRNAU105 


Independent 


CpG island within £ 2Kb 






U104 


CDBox 


17 




30 


28S rRNA C1327 


Independent 


CpG island within £ 2Kb 






UGO 


CDBox 


16 




83 


28S rRNA G4340 


Independent 


CpG island within £ 2Kb 






U3 


CDBox 


1 7 




2 17 


Unknown 


Independent 


CpG island within £ 2Kb 






U3-2 


CDBox 


17 


* 


217 


Unknown 


Independent 


CpG island within £ 2Kb 






U3-2B 


CDBox 


17 




217 


Unknown 


Independent 


CpG island within £ 2Kb 






U3-3 


CDBox 


17 




217 


Unknown 


Independent 


CpG island within £ 2Kb 






U3-4 


CDBox 


17 




217 


Unknown 


Independent 


CpG island within £ 2Kb 






U98b 


HAcaBox 


1 




133 


Unknown 


PPP2R5A 


Host gene with a 5'CpG island 






ACA20 


HAcaBox 


6 




132 


18S rRNA U651 


TCP1 


Host gene with a 5'-CpG island 






ACA29 


HAcaBox 


6 




140 


Unknown 


TCP1 


Host gene with a 5 -CpG island 






U32A 


CDBox 


19 




77 


18S rRNA G132B and 28S rRNA A151 1 


RPL13A 


Host gene with a 5-CpG island 






U33 


CDBox 


19 




83 


18S rRNA U1326 


RPL13A 


Host gene with a 5'-CpG island 






U34 


CDBox 


19 


* 


66 


28S rRNA U2824 


RPL13A 


Host gene with a 5 -CpG island 






U35A 


C DBox 


19 


* 


86 


28S rRNA C4506 


RPL13A 


Host gene with a 5'-CpG island 






SNORD123 


CDBox 


S 


+ 


70 


Unknown 


LOC100505806 


Host gene with a 5 -CpG island 






U70C 


HAcaBox 


9 




136 


18S rRNA U1692 


ASTN2 


Host gene with a 5 -CpG island 






ACA59B 


HAcaBox 


17 


+ 


1S2 


Unknown 


SLC47A1 


Host gene with a 5'-CpG island 






U3BA 


CDBox 


1 


+ 


71 


28S rRNA A185S 


RPS8 


CpG island within £ 2Kb 






U38B 


CDBox 


1 


* 


'3 Ei 


28S rRNA A1853 


RPS8 


CpG island within £ 2Kb 






HBII-52-1 


CDBox 


16 


* 


B2 


serotonin receptor 5HT-2C mRNA? 


StJURF-SI JRNP-UBE3A antisense 


CpG island within £ 2Kb 






HBII-35-3 


CDBox 


15 


+ 


97 


Unknown 


SUURF-SI JRNP-UBE3A antisense 


CpG island within £ 2Kb 






SNORD116 1 


CDBox 


1 




92 


Unknown 


USH2A 


CpG island within £ 2Kb 






ACA66 


HAcaBox 


1 




137 


1SS U3G 


PABPC4 


CpG iolond within £ 2Kb 






USE 


scaRna 


2 


* 


266 


Us snRNA U41 


ATG15L1 


CpG island within £ 2Kb 






HBII-316 


CDBox 


2 


* 


89 


28S rRNA A3846 


WDR43 


CpG island within £ 2Kb 






ACA58 


HAcaBox 


3 




137 


28S rRNA U3823 


MRPL3 


CpG island within £ 2Kb 
















28S rRNA U3741 and 23S rRNA U3743 


jiTP 6 VUE 


CpG island within £ 2Kb 






UG3 


CDBox 


S 




6B 


28S rRNA A4541 


HSPA9 


CpG island within £ 2Kb 






ACA54 


HAcaBox 


11 




123 


28S rRNA U3801 and 23S rRNA U4S39 


NAP1L4 


CpG island within £ 2Kb 






ACA49 


HAcaBox 


12 




13i= 


Unknown 


EP400 


CpG island within £ 2Kb 






ACA2b 


HAcaBox 


12 




137 


23S rRNA U4263 and 28S rRNA U4282 


C12orf41 


CpG island within s. 2Kb 






ACA28 


HAcaBox 


14 




127 


13S rRNA U815 and 13S rRNA U86S 


EIF-3 


CpG island within £ 2Kb 






HBII-82 


CDBox 


16 




M 


28S rRNA G3923 


SF3B3 


CpG island within s. 2Kb 






U105B 


CDBox 


19 




79 


18S rRNA U799 


PPAN 


CpG island within s. 2Kb 






U37 


CDBox 


19 






28S rRNA A3697 


EEF2 


CpG island within £ 2Kb 






HBII-180A 


CDBox 


19 




97 


28S rRNA C3680 


C19orf48 


CpG island within £ 2Kb 






HBII-180B 


CDBox 


19 






23S rRNA C3630 


C19orf48 


CpG island within £ 2Kb 







Green and red rectangles represent unmethylated and methylated CpG islands, respectively. 



microRNAs 14 and T-UCRs 18 undergoing cancer-specific CpG 
island hypermethylation-associated silencing are known to have 
tumor suppressor roles, it is possible that snoRNAs act in a simi- 
lar manner. This additional level of complexity is really true 
for the epigenetic silencing of two of the identified snoRNAs, 
SNORD123, 50 and ACA59B (also known as SNORA59B), 51 
because their target RNAs are unknown. ACA59B resides in an 
intron of the Solute Carrier Family 47 Member 1 (SLC47A1) 
gene, while SNORD123 is a C/D box snoRNA that resides 
within a long ncRNA transcript (LOC100505806) while in 
opposite direction is transcribed from the same CpG island the 
coding gene SEMA5A, adding another level of complexity in 
this case. For U70C (also known as SNORA70 or ACA70), the 
task might be a little easier. U70C was originally cloned from 
a cervical cancer cell line and belongs to the H/ACA box class 
of snoRNAs, having the predicted hairpin-hinge-hairpin-tail 
structure, conserved H/ACA-box motifs, and an association 
with the GAR1 protein. 50 ' 52,53 The snoRNA U70C resides in an 
intron of the Astrotactin 2 (ASTN2) gene in the sense orienta- 
tion and it serves as a guide for the pseudouridylation of selected 
bases of rRNA by forming short duplexes with the 18S rRNA 
U1692, the target for this snoRNA. 50 ' 52 ' 53 A role for 18S rRNA in 



tumorigenesis is starting to emerge, 54 " 56 and our findings provide 
additional information about this role. 

The enormous task of understanding the mechanisms by 
which snoRNA epigenetic silencing contributes to the origin 
and progression of human tumors still lies ahead. In the mean- 
time, our observation that epigenetic inactivation by CpG island 
hypermethylation of a subset of snoRNAs, such as SNORD123, 
U70C and ACA59B, occurs across a wide spectrum of human 
cancer cell lines and primary tumors of diverse cellular and tissue 
origin provides clear support for the concept that major disrup- 
tion of ncRNA programming is a common feature of cancer cells. 

Materials and Methods 

Cell lines, culture conditions and primary study samples. The 

human cancer cell lines examined in this study were obtained from 
the American Type Culture Collection (ATCC). HCT-116 and 
DKO cells were a generous gift from Dr. Bert Vogelstein (Johns 
Hopkins Kimmel Comprehensive Cancer Center). Cell were 
maintained in appropriate media and treated with 1 u,M 5-aza- 
2'-deoxycytidine (Sigma) for 48 h to achieve demethylation. 18 
DNA samples from normal tissues and primary leukemias were 



884 



RNA Biology 



Volume 9 Issue 6 



PPP2R5A (U98b) 

Normal Colon 

I I II I I III 1 1 h 



■ 



inn iiiii ^M iiii i i iiii m i i 3 



TCP1 (ACA20 and ACA29) 

, Normal Colon 



5'. 



i— m — m i i i iii i ii iiiiii hhw- 



+131 -280 



■H-H — I III II I I I II 



BEtEBB BBS 



H 1 



+258 



s HCT116 

H—m HW h 

-397 I 



i hi m i i — in i iiiii 1 1 ii i 1 1 i i h i n i i n i n — i 



g HCT116 

I — HH — H-H- 



— H-t+d- 



■H-H — I II I I I I I II 



H 1 



SNHG5 (U50 and U50B) 

g Normal Colon 3 . 

I Hfl — III III M l IIIIII II I II III I IH-H — I — H — H 1 1 II I I I 1 



RPL13A (U32A, U33, U34 and U35A) 

;> Normal Colon 



I I I I I — I 1 — I — HH II II M il II III I I I I HhH HHH— I 1 1 I I I — I 



+541 -442 I 



HCT116 



S 
I— 
+85 



Ht 


M\\ 




III III II 


IIIIII II 1 II 


-H 


\ 

■ 


41+ 


H 


\ — 


1 — H— 

■ 


++ 


— i 


— h 
















] 








! [ 


1 











HCT116 



+541 -442 I 



1 II I I 

+46 



Figure 2. Bisulfite genomic sequencing of the CpG islands associated with the snoRNAs U98b (host gene PPP2R5A), ACA20/ACA29 (host gene TCP1), 
U50/U50B (host gene SNHG5) and U32A/U33/U34/U35A (host gene RPL13A) in normal colon and the colorectal cancer cell line HCT-116. CpG dinucleo- 
tides are represented as short vertical lines. Eight single clones are represented for each sample. Presence of a methylated or unmethylated cytosine is 
indicated by a black or white square, respectively. Transcription start sites are represented by vertical black arrows. 
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Figure 3. Bisulfite genomic sequencing of the CpG islands associated with the snoRNAs U88, U19-2, U63 and ACA49 in normal colon and the colorectal 
cancer cell line HCT-116. CpG dinucleotides are represented as short vertical lines. Eight single clones are represented for each sample. Presence of a 
methylated or unmethylated cytosine is indicated by a black or white square, respectively. Transcription start sites are represented by vertical black 
arrows. 



obtained at the time of the clinically indicated surgical procedures. 
All patients gave written consent to participate in the study and the 
Ethics Committee of each hospital approved the study protocol. 

DNA methylation analyses. The CpG Island Searcher 
Program 45 was used to determine which snoRNAs were located 
within 2 Kb of a CpG island. DNA methylation status was estab- 
lished by PCR analysis of bisulfite-modified genomic DNA, 
which induces chemical conversion of unmethylated, but not 



methylated, cytosine to uracil. Two procedures were used. First, 
methylation status was established by bisulfite genomic sequenc- 
ing of the corresponding CpG islands. Eight independent clones 
were analyzed. The second analysis used methylation-specific 
PCR with primers specific for either the methylated or modified 
unmethylated DNA. The primers used are described in Table SI. 

Quantification of snoRNAs with real-time PCR. 
Quantitative real-time PCR was performed to quantify the level 
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Figure 4. Bisulfite genomic sequencing of the CpG islands associated with the snoRNAs SNORD123 (host gene is the IncRNA L0C1 00505806), ACA59B 
(host gene SLC47A1) and U70C (host gene ASTN2) in normal colon and the colorectal cancer cell line HCT-116. CpG dinucleotidesare represented as 
short vertical lines. Eight single clones are represented for each sample. Presence of a methylated or unmethylated cytosine is indicated by a black or 
white square, respectively. Transcription start sites are represented by vertical black arrows. 



of snoRNAs, as described previously in reference 18. Briefly, to 
quantify SEMA5A and the snoRNA-host genes 1 (jug of puri- 
fied and DNase-treated (TURBO DNA-free, Ambion) total 
RNA was reverse-transcribed using Thermoscript RT and ran- 
dom primer hexamers. cDNA was amplified by real-time PCR 
using SYBR (Applied Biosystems) green detection. Reverse 
transcription using a custom-designed TaqMan microRNA 
Reverse Transcription Kit (Applied Biosystems) was used to 
quantify the SNORD123 and U70C, providing specificity for 
the mature RNA target. Reactions were performed in an Applied 
Biosystems 7900HT Fast Real-Time PCR system in 384-well 
plates. Expression values of ASTN2, LOC100505806, SLC47A1 
and SEMA5A were normalized to HPRTl and expression values 
of U70C and SNORD123 to RNU19, respectively. Total RNA 
was extracted from three independent experiments and real-time 
PCR reactions were performed in triplicate. The primers used are 
described in Table SI. 

Northern-blot. Fifteen micrograms of total RNA were 
loaded in a 15% denaturating polyacrylamide gel containing 8 
M Urea in 0.5x TBE buffer system. Decade Marker (Ambion) 



was prepared according to manufacturer's protocol, using ["y- 32 P] 
ATP (PerkinElmer) and simultaneously loaded into the gel. Both 
RNA and marker were resolved by gel electrophoresis and trans- 
ferred onto Hybond-N* membrane (Amersham) in 0.5x TBE, 
followed by UV-cross linking (1,200 Jules). Both SNORD123 
and 5.8S rRNA probes were radiolabeled with 25 (JtCie using T4 
kinase (Invitrogen) and purified with Nucaway Spin columns 
(Ambion). The membrane was prehybridized in hybridization 
buffer for 1 h and hybridized overnight in the same solution at 
45°C containing the SNORD123 probe previously heated at 
95°C for 2 min. The membrane was washed at low stringency 
followed by film exposure. The membrane was then hybridized 
with the 5.8S rRNA probe using the same conditions followed by 
quantification using phosphorimager technology. All the probes 
used are described in Table SI. 

Infinium 45 OK DNA methylation array. The DNA meth- 
ylation levels at 10 CpG sites encompassing the SNORD123- 
associated CpG island were determined using the Infinium 
450K DNA methylation microarray, as previously described 
in reference 49. Briefly, DNA was quantified by Quant-iT™ 
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Figure 5. Expression analyses of of the transcripts derived from the 
SNORD123/LOC100505806/SEMA5A CpG island. (A) Quantitative RT-PCR 
of SNORD123, LOC10050580 and SEMA5A showed loss of expression in 
the CpG island hypermethylated HCT-116 and SW48 cells. SNORD123 
ACA59B and U70C are expressed in the unmethylated DLD1 cancer 
cells and in normal colon mucosa. (B) Northern-blot analysis shows the 
absence of the SNORD1 transcript in the hypermethylated HCT-116 and 
SW48 cells and its presence in the unmethylated DLD1 cells. 



PicoGreendsDNA Reagent (Invitrogen) and the integrity was 
analyzed in a 1.3% agarose gel. Bisulfite conversion of 600 ng 
of each sample was done according to the manufacturer's recom- 
mendations for the Illuminalnfinium Assay. Effective bisulfite 
conversion was checked for three controls that were converted 
simultaneously with the samples. Four (jlI of bisulfite-converted 
DNA were used to hybridize on an InfiniumHumanMethylation 
450 BeadChip, following the Illuminalnfinium HD Methylation 
protocol. The chip was analyzed using an IlluminaHiScan SQ 
fluorescent scanner. The intensities of the images were measured 
using GenomeStudio (2010.3) Methylation module (1.8.5) soft- 
ware. The methylation score of each CpG is represented as a (3 
value. 
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Figure 7. SNORD123, U70C and ACA59B CpG island methylation in human cancer cell lines. DNA methylation microarray analyses of five, two and four 
CpG sites corresponding to the bisulfite genomic sequenced CpG islands for SNORD123, U70C and ACA59B in the NCI60 panel of seven different tumor 
types. Each square represents a single CpG: green square, unmethylated CpG; red square, methylated CpG. The three snoRNA-associated CpG islands 
were methylated in the originally studied HCT-116 cells, but hypermethylation events were also observed in other classes of malignancies, such as 
leukemias. 
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Figure 8. SNORD123, U70C and ACA59B CpG island methylation in primary acute lymphoblastoidleukemias. DNA methylation microarray analyses 
of five, two and four CpG sites corresponding to the bisulfite genomic sequenced CpG islands for SNORD123, U70C and ACA59B in leukemia patients 
demonstrated hypermethylation. Each square represents a single CpG: green square, unmethylated CpG; red square, methylated CpG. The B-cell or 
T-cell phenotype of the studied ALLs is indicated. 
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